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Abstract 

An approach to the electrostatic control of 90° magnetization rotation in the hybrid structures 
composed of topological insulators (TIs) and adjacent ferromagnetic insulators (FMI) is proposed 
and studied. The concept is based on TI electron energy variation with in-plane to put-of plane 
FMI magnetization turn. The calculations explicitly expose the effect of free energy variability in 
the form of the electrically controlled uniaxial magnetic anisotropy, which depends on proximate 
exchange interaction and TI surface electron density. Combining with inherent anisotropy, the 
magnetization rotation from in-plane to out-of-plane direction is shown to be realizable for 1.7 
~ 2.7 ns under the electrical variation of TI chemical potential in the range ± 100 meV around 
Dirac point. When bias is withdrawn a small signal current can target the out-of-plane magneti- 
zation instable state to the desirable direction of in-plane easy axis, thus the structure can lay the 
foundation for low energy nonvolatile memory prototype. 

PACS numbers: 75.70.Cn, 75.75.Jn, 73.63.Rt, 85.70.Ay 
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Because of the absence of the natural magnetic charges and thus the leakage conditions, 
magneto-electronic automata presents a great advantage of the magnetic devices compared 
with their electrical counterparts in terms of the non- volatility. Up to date the magnetization 
driving force in such devices is realized via a magnetic field of the channel current or current- 
induced angular momentum transfer in magnetic tunnel junction.- The current-assisted data 
processing in both techniques consumes much energy, the reduction of which is one of the 
major concerns for present industrial developments and academic researches. 

Electric field control in addition to the current induced effects would has been pursued 
with multiferroic materials. The tension variation in some FM materials can rotate their easy 
axis 90° jointly with magnetization vector M.-~- In FM/piezoelectric hybrid structures such 
effect evokes the new paradigm of magnetization rotation in the electrical field. However it 
should be noted that the desirable piezo-effect in such structures requires large gate electric 
field, which indicates the need of searching for much stronger magnetoelectric effects in 
modern magneto-electronics. 

Another approach relies on the effective magnetic fields H ex mediated by the proximate 
exchange interactions between surface magnetic ions and itinerant electrons (or holes) in 
semiconductors,- graphene,- - - or topological insulators.— 1 ^ A distinctive feature of the H ex 
consists in its dependence on electron density at the interface of the electric channel and 
magnetic layers^ 1 ^ rather than on an electric current. This property of the H ex combines 
the magnetic and electric responses inside the new "meta-materials" in a strongly correlated 
manner that stems from the electrostatic and quantum mechanical nature of the exchange 
interaction. 

Recent studies have demonstrated the electrical control efficiency of the proximate ex- 
change interaction (PEI) between a magnetic layer and graphene for different spintronic 
applications.- 1 ^ 1 ^— As an example, the conductivity of the graphene sandwiched between 
ferromagnetic insulators (FMIs) reveals the giant magnetic resistance.— 1 ^ On the other 
hand, the magnetization of FMI can be turned and even switched by electrical variation of 
the carriers density in adjacent graphene.— 

In present letter we explore the new prospects granted by topological insulators (TIs) 
and composite structures (Fig. 1) based on TIs. Specific spin-momentum coupling in 3D 
TI's surface Dirac fermions can be naturally incorporated into magneto-electric effects of 
TI/FMI composite structures.— As a result, a lot of unusual physical phenomena have been 
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predicted in TIs affected by adjacent FMIs.— — Here we focus on converse TI effect on 
adjacent FMI. Namely, the PEI of TI surface electrons is shown to appear in the form of FMI 
effective magnetic anisotropy electrically driven without any structure tension. Combining 
with inherent anisotropy and anisotropy of the sample shape, the in-plane to out-of-plane 
orientation transition can be electrically stimulated by variance of Tl-mediated anisotropy. 
This effect is quantitatively demonstrated for uniaxial nanomagnet deposed on the surface 
of TI (e.g. Bi2Ses). It is remarkable that the system with perpendicular M appears in 
neutral position with respect to M returning to either of the in-plane states along easy axis. 
This instability can be used for directional targeting of M by a current induced torque.— 
Small signal current is sufficient to tilt the magnetization to the desired direction and the 
following switch will be driven by the shape and intrinsic anisotropy terms. The possible 
design of nonvolatile memory and logic devices^ can be realized on basis of these magnetic 
units. 

Qualitatively the easy axis rotation can be thought as an effect of the electron energy 
reducing under exchange interaction with proximate FMI. When the magnetization is in 
the in-plane direction, the PEI results in a shift of the Dirac cones in the momentum space 
whereas an energy gap can be generated and the energy dispersion becomes non-linear 
for the out-of-plane orientation with M z ^ 0.— From the thermodynamic point of view, 
the case with a band gap (M z ^ 0) is expected to be energetically more favorable with 
lower electronic free energy than that with the shifted bands (M z = 0) following a similar 
consideration demonstrated in bilayer graphene.— A departure of chemical potential \x from 
the energy level of the Dirac point Ed suppresses this tendency because of energy increase 
of the conduction electrons partially compensates valence band energy decrease while gap 
opening for ji > Ed and fewer active electrons for \i < ejy- 

To proceed the quantitative analysis let us consider a thermodynamical potential of TI 
electrons interacting with FMI, 

= -T J> (l + exp »- £D - £b A (1) 

6,k ^ ' 

where the sum involves electronic bands b and wave vector k, T is the temperature in energy 
units. We are interested in alteration of $ e (/f) with the Fermi level varying near the Dirac 
point of the surface states. Thus the rest of electronic states in Eq. ([[]) are irrelevant and 
summation in $ e (/i) can be safely restricted by the spectrum of surface Dirac fermions with 
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2D effective Hamiltonian^ 



H = hvp [a x k] z + Dk 2 + Ga (2) 

where vp is the electron Fermi velocity near Dirac point, a the Pauli matrix of electron spin, 
D is the material parameter for quadratic term. The remaining term in Eq. (J2J) describes 
the energy of an electron spin in the exchange effective fields (in units of energy) G = Gm 
of the proximate FMI, m = M/ |M| ensures the collinearity of the effective field and FMI 
magnetization. Introducing the angles 9 between z and m and ip between k and planes (m, 
z) (Fig. 1) the energy spectrum of Hamiltonian (|2J) can be expressed as 

£b,k = D p 2 + b Vp 2 + G 2 + 2Gp sin 9 sin <p, (3) 

where p = fa^k is the electron surface momentum in energy unit, Dq = D/h 2 v F and b = ±1 
marks the conduction and valence bands. Equation (jHJ) clearly illustrates the band structure 
reconstruction lowering the energies of electrons near the top of valence band when M turns 
from in-plane direction {9 — tt/2) to vertical position (9 = 0). The valence band electrons 
(6 = —1) contribute to this effect while conduction band electrons (b = 1) may reduce 
it. Thus the proximity with TI plays the role of a mediator enhancing the uniaxial FMI 
anisotropy along the normal z. 

Following Eqs. (j3J) and ([I]), $ e (A i ) is an even function of G x = Gsin9 that determines the 
dependence $ e (/i) = —K e s(fi) sin 2 9 and proportionality K e g(fj,) ~ G 2 in low-order expansion 
on G. Similar relation remains for A$ e (/i) = $ e (££>) — $ e (/i), which reflects the variable 
part of the TI mediated anisotropy energy 

A$(/i) = AAK cS (/2) sin 2 9. (4) 

Here A is the area of FMI/TI interface and 

AK cS (fi) = f(fi,T)D G 2 , (5) 

where Do is introduced for convenience sake. Computation of Eq. (CQ) with Eq. fl3]) confirms 
the validity of the definition in Eqs. (]1J) and (jSJ). Function /(/x, T) can be numerically found in 
terms of Eqs. (CD) and ([3]); Fig. 2 depicts the results of calculations. The results indicates the 
amplitude of anisotropy changes with surprisingly minor temperature effect but significantly 
asymmetry in chemical potential variation over conduction band (/i — Ed > 0) and valence 
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band (fi — Eo < 0). Consequently, Eq. dSJ) predicts the variation range of anisotropy as large 
as 0.5 meV/nm 2 provided the strength of exchange interaction is G = 40 meV.— In such a 
case a contact area section of tens nanometers may supply the anisotropy energy change of 
eV-scale that guarantees the nonvolatile magnetic bi-stability at room temperature. 

Whether the variation of AK e g(fi) is capable to cant over M depends on symmetry and 
magnitude of total magnetic energy F of the nanomagnet. It is convenient to approximate 
the actual shape of the FMI (e.g. rectangular parallelepiped) by triaxial ellipsoid with 
homogeneous magnetization over the whole volume V . At such approach the free energy 
can be simplified to normal form 



where first term describes the energy of demagnetizing field, and last one represents the 
magnetic anisotropy, which includes independent on \x part $ e (££>) of PEI and tension- 
mediated term with uniaxial anisotropy, Nij the components of demagnetizing tensor. In 
the principal coordinates i,j = x,y,z, this tensor becomes diagonal, = Nidij, with 
principal values 



where £\ iVj = 1 and Qi = q + a| with lengths a x ,a y ,a z of the ellipsoid principal axes. If 
they coincide with crystalline axes, the free energy per unit volume can be brought into 
canonical form 



with u an (m) = U an (m) /V . In the case of planar FMI/TI structure (i.e. N z » N x ,N y ) the 
demagnetizing field tends to establish the hard axis along normal z to interface. So, in soft 
magnetic materials like Yttrium Iron Garnet, the proper choice of FMI sizes can ensure the 
electrical switch between in-plane to out of plane M direction. This condition, however, is 
not enough to guarantee the nonvolatile bistability because there is a path that connects 
the states M x and —M x through the saddle points ±M y . Apparently demagnetizing field 
can also establish the in-plane anisotropy in the structures with shorter FMI layer length in 
y-direction, but desirable retention time with saddle point height ~ 2irM 2 (N y — N X )V = 1 
eV can be reached at least in sub-micron device sizes (200 -300 nm) as N y is still relatively 
small compared with N z . 




(6) 




(7) 




(8) 
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An example of the alternative approach assumes the FMI with intrinsic uniaxial crys- 
talline anisotropy in form of hard y-axis. In such a case the anisotropy term is 

u an {m) = K im 2 y + [K u + -AKM}(1 - m 2 z ), (9) 

where K\ > 0, K u accumulates the effects of crystal lattice distortion, which, in turn, vary 
with crystal doping, the strain at interfaces with TI and spacer layer attached on opposite 
sides, and PEI accounted for $ e (££>). Thus, this parameter depends on particular device 
implementation. 

The total magnetic energy expressed in spherical coordinate of m (i.e., m x = sin 9 cos <p; m y = 
sin^suK/?; m z = cos 6, see Fig. 1) is shown in Fig. 3(a) and 3(b) for two chemical potentials. 
The graph clearly illustrates the alternation of energy valleys and energy hills located along 
in-plane directions, 9 = n /2 and (p = or n, [Fig. 3(a)] and in the normal directions, 9 = or 
Ti, [Fig. 3(b)], as chemical potential changes. A small perturbation of unstable states at tops 
of energy hills may guide the following relaxation into desirable valley. For such perturbation 
we suggest to evoke a proximity induced effective magnetic field B cff = [J x z] G/evpL z M 
(e is the electron charge, L z is the FMI thickness) which accompany the spin polarized TI 
surface current J.— — This field can be used to choose the path to the desired valley by the 
current direction. Figure 3(c) and (d) illustrates the distortion of energy relief under small 
electric current. All together Fig. 3(c)-(d) exhibit the possibility of the electric field control 
of magnetization turn and even switch provided appropriate handling by bias voltage and 
signal current. 

The dynamics of magnetization turn determines the switching speed. A common ap- 
proach to this problem is based on the solution m = m(t) of the Landau-Lifshitz-Gilbert 
equation, which, in the case of mono-domain approximation with homogeneous magnetiza- 
tion, takes the form: 

dm dm 

— = - 7 m x H e// + am x — . (10) 

where 7 is the gyromagnetic ratio, a the Gilbert damping factor and the effective field 
H c fj = (nQMoV)~ 1 dF / dm stems from the total free energies F [Eq. (JBJ)]. Note that H cff 
crucially depends on the magnetization.— For example, if the initial magnetization m(0) 
perfectly aligned with x axis and ignoring fluctuation, the H e ff and m are collinear and the 
system comes into frustrating (metastable) state with respect to subsequent relaxation to 
the z or — z direction. However any thermal dispersion of m(0) may work as a start point 
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for magnetization switch. To show this effect we consider evolution m(£) starting from 9 and 
ip that deviate from idealistic case 9 = 90° and <fi = 0° on only 0.1°. The resulting switch 
is shown in Fig. 4. The switch is gradually accelerated till t = 1 ns, after which the H c g 
gathers strength and steep switch occurs with oscillation. The actual switch time is about 
1.7 ns marked by m z changing from 0.1 (point B) to 0.9 (point A), while it may take around 
1 ns for the switch to be initiated. 

To see the dispersion of switch time, we initializes the random array of magnetization 
states according to Boltzmann distribution at room temperature. The dispersion of the shift 
is marked by the square point A (when m z ^ 0.9) in Fig.4 and the range is indicated by 
the corresponding arrow. This means the total 90° switch time varies from 1.7 ns to 2.7 
ns because of thermal fluctuation. One possible solution to the undetermined switch time 
is to use a TI surface current to initiate the switch by current induced torque. Thus 180° 
switches can be achieved. 

In conclusion, a magneto-electric effect in form of electrically controlled magnetic 
anisotropy has been presented for the structure taking advantage of the proximity ex- 
change interaction between a TI and a FMI layer. The effect is based on exchange energy 
correlation with weak temperature dependence and it is more respondent for the change 
in valence band than that in the conduction band. This interaction effectively induces a 
change of the perpendicular anisotropy energy of the FMI layer which results in a thermally 
initiated magnetization rotation from in-plane to out-of-plane. The initiating time can be 
as long as 1.4 ns if well aligned and the switch time is about 1.7 ns. A small signal current 
can be applied to both initiating the switch and to determine the magnetization after with- 
drawing the bias by tilting the magnetization to the desired direction at the beginning of 
the relaxation. 

This work was supported, in part, by the SRC Focus Center on Functional Engineered 
Nano Architectonics (FENA) and the US Army Research Office. 
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FIG. 1. (Color online) Schematic illustration of the structure composed of ferromagnetic insulator 
layer (FMI) deposed on topological insulator (TI). The controlling gate is not shown. The x,y,z 
are the reference frame used in the paper. The easy axis (two-sided arrow) is along x-direction 
without gate voltage. The 6 and tp are the spherical coordinates of the off-balance magnetization 
vector M. A signal current can be applied to the TI surface to tilt the magnetization. 
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FIG. 2. (Color online) Parameter /(/x,T) reflecting the TI mediated anisotropy as a function of 
chemical potential \i at three temperatures 300 K (curve 1), 200 K (curve 2) and 77 K (curve 3). 
Calculations use the parameters of Bi2Se3 for Hamiltonian : vp = 6- 10 cm/s; D = 13.0 eV- A 2 . 
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FIG. 3. (Color online). Energy landscape in polar coordinates 6 (radial variable) and ip (azimuthal 
variable, Fig. 1) of magnetization direction. The equal-energy curves are calculated for sizes of FMI 
planar layer 50 x 50x2 nm , M = 4.7T x 160 G, K x = K u = 0.05 meV/nm 3 (a) Chemical potential 
and signal current are fi — Ed = —90 meV and J = 0. Two symmetrical valleys correspond to M 
orientations along easy axis x. (b) Easy axis is turned to z-direction with /i — Eo = and J = 0. 
(c) Signal current shifts the valleys and tilts magnetization to —x direction(/i — ££> = meV and 
J = 1 /iA/nm) that drives the M relaxation to the reversal direction pointed at figure (d) with 
(fi — Ed = —90 meV and J = 1 /iA/nm). 
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FIG. 4. (Color online) Anisotropy driven switching trajectory. The initial magnetization is set to 
9 = 89.9°, 4> = 0.1° to enable the initiation of the switch. Square B (m z = 0.1) and square A 
(m z = 0.9) mark the starting and ending of the switch respectively. The switching time is around 
1.7 ns while the initiating time can be as long as 1.4 ns due to thermal fluctuation at 300 K.The 
finishing point A thus varies in the range indicated by the double arrow. 
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